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Abstract 

A shear-free ray congruence on Minkowski space is a 3-parametcr fam- 
ily of null geodesies along which Lie transport of a complementary 2- 
dimensional spacelike subspace (called the screen space) is conformal. 
Such congruences are defined by complex analytic surfaces in the associ- 
ated twistor space CP^ and are the basis of the construction of massless 
fields. On a more general space-time, it is unclear how to couple the 
massless field with the gravitational field. In this article we do this by 
considering the following Cauchy-type problem: given a Riemannian 3- 
manifold {M^,go) endowed with a unit vector field Uq that is tangent 
to a conformal foliation, we require that the pair extend to a space-time 
(M., Q) endowed with a spacelike unit vector field Ut in such a way that 
Ut simultaneously generates null geodesies and is tangent to a conformal 
foliation on spacelike slices t — const. 
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1 Introduction 

Twistor theory, as developed by R. Penrose and others over the last forty 
years, provides a fundamentally different perspective to the quantisation of 
physical fields, by considering the twistor space, the space of null geodesies, 
as the more basic object from which space-time points should be derived. 
The program has had much success as well as some difficulties, particularly 
in its inability to adapt in any natural way, to general curved space-times. 
The article [12] provides background and motivation to the subject. 

One of the appealing aspects of the twistor approach, is the way complex 
analytic geometry comes into play. Specifically, by what is essentially the 
Klein-Pliicker correspondance, a light ray in Minkowsi space corresponds 
to a point in a 5-dimensional CR-submanifold of CP^. To make the cor- 
respondance bijective, one needs to conformally compactify Minkowski space 

4 

by adding a light-cone at infinity to obtain the manifold M diffeomorphic 
to 5^ X The picture can be unified by introducing the flag manifold F12 
of pairs {i, 11) consisting of resp. 1- and 2-dimensional subspaces of with 
i CiL and considering the double fibration: 

F12 

/ \ 
CP3 G2(C^) 

u u 
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where G2{C^) is the Grassmannian of complex 2-dimensional subspaces of 
and where the left projection is given by {£, H) i— > ^ and the right by 

{i, n) ^ n. 

In general, a complex analytic surface S C CP'^ which intersects 
defines a 3- (real) parameter family of geodesies in a region A C (called 
a congruence) which is without shear. That is, if W represents the tangent 
vector field to the congruence, then at a point x £ A, the metric complement 
W-^ is 3 dimensional and contains W itself; if we take a 2-dimensional 
spacelike complement S in W-^, then for the congruence to be shear-free, 
we require Lie transport of vectors in S along W to be conformal. Such 
congruences can be used to generate solutions of the massless field equations. 
This arises by contour integration in the twistor space and is the basis of the 
Penrose transform, which is an integral transform from sheaf cohomology in 
the twistor space into the space of massless fields, see [8, 16] for details. 

If we consider Minkowski space endowed with standard coordinates 
(t = , , x'^ , x"^) , then in [3] it was shown how an SFR corresponds to 
a semi-conformal mapping (see Section 3 for the definition) on a space-like 
slice Rj := {{t, x^ , x^ , x^) G : t constant}. The argument given in [3] is 
rather indirect and depends on the analyticity of the SFR; it is essentially 
as follows. An SFR defines a complex analytic surface S in CP^. If we 
now regard CP^ as the (Hopf) bundle whose fibre at each point x £ S'^ is 
the set of almost Hermitian structures at x, then S determines an integrable 
complex structure J on a subset V oi S'^. If we consider as being made up 
of two 4-balls joined along an equatorial 3-sphere C 5^, then on endowing 
one of these balls with the Poincare metric of constant curvature —1, by a 
theorem of J. C. Wood [18], the complex structure J determines a harmonic 
morphism on V (cf. [1, 4]) which extends to a semi-conformal mapping on 
VCidB^ [3]. Provided one is careful with conventions, dB^ = corresponds 
to the space- like slice Rg. In Section 3, we give a direct argument that an 
SFR corresponds to an evolving family of complex-valued semi-conformal 
mappings on space- like slices R^. 

This point of view is readily adapted to a more general construction: we 
propose to replace the SFR with an evolving family of unit vector fields Ut 
on a 3-manifold which is coupled to a family of Riemannian metrics gt 
on M^. We require that the pair (gt, Ut) generate a space-time metric of the 
form Q = —fdt"^ + gt, for some function / G C°°{I x M^) (/ C R), in such a 
way that Ut determines a null vector field dt + fU which is both geodesic and 
shear-free. Here, we have to be careful about what we mean by shear-free, 
for there are two natural definitions which in general only coincide when 
Q is the fiat metric. For the first of these, we consider a ray congruence 
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along which Lie transport of vectors in a complementary spacelike 2-space 
is conformal, i.e. neighbouring rays are without shear; this is the usual 
notion of a shear-free ray congruence (SFR) as discussed above in the case 
of Minkowski space. 

An alternative notion is to require that for each t, the vector field Ut 
be tangent to a conformal foliation (or shear-free foliation) on {M^,gt). In 
order to formulate a well-posed problem, we still require that Ut generate 
null geodesies and will refer to the corresponding congruence as one that 
generates conformal foliations (CFGR: "conformal foliation generating ray 
congruence"). It should be remarked that if a vector field U is tangent to 
a conformal foliation on a domain A of a Riemannian 3-manifold {M^,g), 
then we can locally integrate it to a semi-conformal map (p : V C A ^ C, 
defined up to post-composition with a conformal tansformation. Thus the 
vector field Ut and the map ipt are locally interchangeable. Furthermore, 
a semi-conformal map is absolutely minimizing for a natural functional, 
and we will see in Section 6 that the family {(ft} of semi-conformal maps 
associated to the family {Ut} satisfy a differential equation somewhat akin 
to the Schrodinger equation. 

Our construction is analogous to the Cauchy problem, usually formulated 
to generate solutions to the Einstein equations, see [5] for an overview. 
However, rather than produce an Einstein manifold, we wish to generate a 
space-time that supports a CFGR (or an SFR) for a given initial condition 
(M^, go, Uq), where Uq is tangent to a conformal foliation on (Af^, (70). Thus, 
instead of defining a field with respect to a fixed background metric, we 
suppose the field ft and the metric gt are coupled so as to produce a CFGR. 
We are led to consider the evolution: 



with the above initial conditions. Here, = g{U, ) is the 1-form dual to 
U and T is a t-dependent symmetric covariant 2 tensor on A!^ satisfying 
the property T{U, U) = at each point. Our key objective is to choose / 
and T in such a way that a solution produces a CFGR on the space-time 
(Al^ = / X M^,Q = — /^di^ + gt). In Theorem 5.2, we give condition when 
this is the case. 

We study in particular, fixed points of the evolution, which occur when 
go is conformally equivalent to a metric with respect to which the vector 
field Uq is tangent to a conformal foliation by geodesies; the case of con- 
stant curvature, when we can generate explicit solutions; the case when the 





fVfjU + gradgf 
26 (d/lTAf) + 2T 
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distibution on M"^ orthogonal to Uq is integrable-remaikahly, this property 
is preserved under the flow and the metric evolves according to a modi- 
fied Ricci flow (Theorem 5.6). We also study a similar coupled evolution 
with given initial data, in order to generate a space-time supporting an SFR 
(Theorem 4.7). In the final section we discuss the functional representation 
of our construction, its relation to physical fields and conformal invariance. 



2 Shear-free ray congruences on Minkowski space-time 

In order to describe light rays in Minkowski space, we adopt the spinor no- 
tation of Penrose [10]. Thus, to a (real) vector (x") = (t = , , x"^ , x^) G 
M^, we associate the Hermitian matrix: 

AA'\ ^ ( t + x^ + ix^ 



In particular, detX = t^ — {x'^Y ~ {x^)'^ ~ {x'^Y' is the Minkowski norm so 
that det X = if and only if x^"^' = ^^^^ , where to the spinor G we 

—A' —0' — —1' — 

associate its conjugate ^ whose components are given by ^ = i^,^ = £}■ 
Indices are raised and lowered using the skew forms: 

{sab) = ie^^) = isA'B') = (e^'^') = ( ; ) 
and as usual we sum over repeated indices. 
The partial derivatives VaA' are defined by 

ryj ^ _ f do + di 82- ids 

where we write dj = d/dx^ , so that Vaa'X^^' = ^a^A' where 6^ and 6^', 
are the Kronecker delta symbols. 

A shear-free ray congruence (SFR) is a congruence of null geodesies such 
that parallel transport of a 2-dimensional complementary spacelike subspace 
(the screen space) along each geodesic of the congruence is conformal (i.e. 
is without shear). Such a congruence is described by a spinor field /x^(x), 
(x G M^) satisfying the equation 

(2) /U^/U^V^A'/iB = . 

Explicitly, at each point x G M^, the vector [/i'^(x), /i^(x)] G CP^ ~ 5^ 
defines the direction of the light ray passing though point in the 

celestial sphere 5^. 
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The equations for an SFR are invariant by a rescaling ^^{x) X(x)fi^{x) 
for some scalar function A, as is apparent when we write (2) in the equivalent 
form: 

^ ' 1 nVov^^ + Vll'^i = 

where /x = fjP / fJ^^- 

Let C CP^ be defined by the homogeneous equation: 

(4) + + ^v""' + W = 

in terms of homogeneous coordinates {(.AtV"^') ^ C^. The correspondance 
between a point [(,a, ] ^ iv^ 7^ 0) ^'^cl points of are given by the 
incidence relation 

(5) Ca = ixAA'V"^' 

This is a matter of linear algebra. In general solutions to (5) will be complex. 
A solution xaa' is real if and only if the corresponding matrix {xaa') is 
Hermitian. Under this assumption, on taking the complex conjugate of 
(5), we obtain four equation in the four unknowns xaa', which we may 
abbreviate by Ax = b, where the coefficients of A are given by rj^ and their 
conjugates and b is given by ^a and their conjugates. Since the determinant 
of A vanishes, there will be either no solutions or infinitely many solutions. 
By Cramer's rule, the latter case occurs if and only if the determinant of 
the matrix obtained from A by replacing the last column by b, vanishes. It 
is straighforward to check that this is precisely the condition (4). Note that 
if XAA' is a real solution to (5), then so is xaa' +^VA'nA'^ for all A G R. This 
is the light ray passing through x^^ with direction r]^ . To complete the 
correspondence, one is required to compactify by adding a light cone at 

4 1 '3 

infinity to obtain the space M diffeomorphic to S x S . Then the points 
[^A,0] G (a copy of CP^) are mapped to null rays generating this cone, 
see [10] for details. 

We can obtain another way. The Hopf fibration vr : CP^ — > is the 
map given by 

7r([zi, Z2, Z3, Z4]) = [Zi + Z2j, Z3 + Zii] G MP^ 

where for clarity, we now use homogeneous coordinates [zi, Z2, z^, Z4] for 
points of CP^ and where MP^ is the quaternionic projective space. On 
identifying MP^ with S"^ and letting be the equatorial 3-sphere given by 
Re [zi + Z2i, zs + z^z]] = 0, we see that vr([zi, Z2, Z3, Z4]) G if and only if 

ZlZ^ + Z2'ZI + + z^za = , 
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which is once more the condition (5). Thus = tt-^{S^). It IS now 
possible to identify 7r~^{S'^) with the unit tangent bundle T^S^ to S^. Then 
the intersection S H J\f^ of a complex analytic surface S in CP^ with J\f^ 
determines a unit vector field on a domain of which is tangent to a 
conformal foliation (see Section 3 for the definition). The basic object which 
defines an SFR is a complex analytic surface in CP^, as can be seen as 
follows. 

Let S C CP^ be a regular complex analytic surface locally parametrised 
in the form {z,w) i-^ w),'!]^ {z,w)], where and rj^ are holomorphic 

in {zjw). Then, provided rj^' ^ 0, up to a relabelling of the canonical 
coordinates of CP^, we can find a reparametrisation of the form {z,w) t-^ 
[a{z,w),T{z,w), (3{z), 1], where now a and r are meromorphic in {z,w) and 
/3 is a meromorphic function of z only. Now the incidence relation (5) takes 
the form 



where we write u = t + x^,v = t — x^,q = x'^ + ix^. On eliminating the 
variable w between these two equations, we obtain a relation of the form: 



where ip is a holomorphic function of three complex variables. Then it is 
quickly checked that any solution z = z{x) of (6) determines a solution 
/i = —l/z of (3) and so defines an SFR. Conversely, any SFR is locally 
given this way; this is know as the Kerr Theorem, see for example [9, 14]. 
Solutions of the zero rest-mass field equations are now given by contour 
integrals of functions ^ given as above. 

Explicitly, the field equations for a zero rest-mass particle of helicity —n/2 
take the form 



where the spinor field ^ab...l is symmetric in its n indices. Special solutions 
are given by contour integrals of the form [11]: 



where / is a holomorphic function of three complex variables. Then we have 
the recurrence relations: 



iV2T 



(5u + q 



(6) 



ip{z, f3{z)u + q, p{z)q + v) = 



V^^Vas...l = 0. 





r = 0,1,2,... 
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For a non-negative integer n, let i^ab—k (with 2n spinor indices) be defined 
by 

V^O = V'OO-OO, V'l = VOO-Ol, ^2 = VOO-ll, • • • , y^AB-K = '^(AB-K) 

Then the recurrence relations above are equivalent to the spinor field equa- 
tion 

These solutions to the zero rest-mass field equations are called algebraically 
special. In their most general form, solutions are obtained via the Penrose 
transform; an integral transform defined on 1-cohomology with coefficients 
in the sheaf of twistor functions homogeneous of degree — n — 2. In general 
the contour in the above integral will surround a pole of / and we obtain 
an SFR via the equation := 1// = 0. 

On writing t for the time coordinate and setting dt = d/dt,di = 
d/dx^,dq = ^{82 — i^a), % = ^{82 + i^s), the equations (3) for a SFR have 
an equivalent expression: 

f (i) fJ-din - fJ^dgiJ, + dgH =0 
\ (u) fldtfl + p?dqn + dq^l = 

These are of course equations expressed in terms of local coordinates and 
there is some arbitrariness in their choice. We first of all show how these 
two equations have an invariant expression, the first in terms of conformal 
foliations on space-like slices Rf := {{t,x^ ^x"^ ,x^) G : t = const} and 
the second as a time evolution of such a foliation. Foliations of this type 
integrate to give locally defined semi-conformal maps, which will be the basis 
of our adaptation of the above constructions to more general space-times. 

3 Reformulation of the SFR equations 

A Lipschitz map ip : {M"^,g) — > {N'^,h) between Riemannian manifolds is 
said to be semi-conformal if, at each point x £ M where f is differentiable 
(dense by Radmacher's Theorem), the derivative d^px ■ TxM — > T^f^^)^ is 
either the zero map or is conformal and surjective on the complement of 
keidcpx (called the horizontal or complementary distribution). Thus, there 
exists a number A(x) (defined almost everywhere), called the dilation, such 
that X{x)'^g{X,Y) = ip*h{X,Y), for all X,Y e (kerdipx)-^. If p> is of class 
C^, then we have a useful characterisation in local coordinates, given by 

where (a;*),(y") are coordinates on M,N, respectively and <pf = i9(y" o 
ip) I dx^ . The fibres of a smooth submersive semi-conformal map determine 
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a conformal foliation, see [15] and conversely, with respect to a local foliated 
chart, we may put a conformal structure on the leaf space with respect to 
which the projection is a semi-conformal map. We then have the identity: 

(Cug) {X,Y) = -2U{lnX)g{X,Y), 

for U tangent and X, Y orthogonal to the foliation. This latter equation can 
be taken to be the characterisation of a conformal foliation. 

The fundamental equation of a semi-conformal submersion (p : M™" 
relates the tension field t((^) := TrVdy? with the gradient of the dilation 
and the mean curvature of the fibres [4]: 

(8) = ~ 2)d99(grad In A) — (m — n)dip{i>) . 

In particular, if n = 2, then a semi-conformal map is harmonic (t((/3) = 0) 
if and only if the fibres are minimal. Note also that any conformal foliation 
only determines a corresponding (locally defined) semi-conformal map up to 
post-composition with a comformal mapping. That is, if ip : — > A^" is 
semi-conformal and ip : N"^ is a submersive conformal mapping, then 

the composition ip o ip : M"^ — > is also semi-conformal with the same 
regular fibre components as ip. 

In order to understand geometrically equations (7), we first prove the 
following lemma. Given a point x of R^, for some t and a unit vector 
U{x) £ TxRf, we define the (future pointing) null geodesic determined by 
U{x) to be the null ray passing through x generated by the null vector 
dt + U{x). 

Lemma 3.1. Let Uq be a unit vector field on a domain o/Rq, which defines 
a vector field Ut on a domain o/R^ by parallel transport along null geodesies 
determined by Uq. Then the evolution of U satisfies 

dU r3 

w = -^-'^- 

Conversely, any solution to this equation with given initial data U\t=o = Uq 
on Rg arises this way. 

Proof. By construction Ut(y) = Uo{x) with y = x + tUo{x). At a point 
X G Rq, we have 

dU 
'dt 



^^Ut{x) - Uoix) 



t=o t^o t 
Now 

Uo{x) = Ut{x + tUo{x)) 

= Ut{x) + tdUt{x)iUo{x)) + 0{t^) 
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so 

dU Uo{x) - tdUt{x){Uo{x)) - Uo{x) +0{t^) 

— — = iim 

dt t^o t 

= -dUo{x)iUo{x)) = -Vu,i^)Uo . 

The converse follows from the uniqueness of a solution to a first order partial 
differential equation, given initial data on a hypersurface. □ 

Recall that a complex- valued function on a domain Minkowski space de- 
fines an SFR if and only if it satisfies equation (3). 

Proposition 3.2. Let fi define a shear-free ray congruence on a domain 
W CM'^ and let 

be the associated unit direction field. Let g = (dx^)^ + (dx^)^ + (dx^)^ he the 
restriction of the canonical metric on to each spacelike slice : t = 
const. Then on W : 




for some function a : ^ R, where X,Y ^ ^R-t o,'^^ vectors orthogonal 
to U and where denotes the Levi-Civita connection on (Rf,g). In 

particular, U is tangent to a conformal foliation on W n Rf , for each t. 

Conversely, any solution U of (9)(i) and (ii) on a domain W C 
determines an SFR on W. Furthermore, if (ii) holds on and (i) holds 
on an initial hypersurface R^y, for some t^, then (i) holds for all t, where 
it is understood that the domain of U may need to be restricted to avoid 
singular points of the associated foliation. Thus conformality of the foliation 
associated to the unit vector field U is preserved under the flow (9)(ii). 

Proof. First note that equation (9)(i) is invariant under rescaling U i— > fU 
for any real valued function /, so we work rather with the vector field 

V = - 1, 2^) = - i)di + 2fidg + 

where we recall that we write = x^,q = x"^ + ix^ . Now observe that (9)(i) 
is equivalent to the equation 

(Cvg) {X-iY,X-iY)=0 

where X,Y £ TRf satisfy ||X|| = ||y||, g{X,Y) = g{X,V) = g{Y,V) = 0, 
and we extend tensors to complex vectors by complex linearity. In what fol- 
lows we use the bracket ( • , • ) where convenient to express the inner product 
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A complex vector with the above properties is given by 

Then, on noting that {dq, dq) = {dq, dq) = and that {dq, dq) = |, we have 

^ (Cvg) {X -iY,X- iV) = giVx^ivV, X-iY)= g{V, Vx~iy{X - iY)) 

= (d^P - 1)^1 + 2^5g + 27^%, V(^a^_^2a^+a_)(^f9i - n'^dq + dq)) 

= -(1 + l/ip) {fldlfl - l?dq^l + dqfl) . 

The equivalence between (7)(i) and (9)(i) now follows. 

We perform a similar calculation in order to determine the time evolution. 
Suppose that (9)(i) and (7)(ii) hold. Then 

V^V = -l)dui + 2f,dqfM + 2jldqfl}i-pdl+2dq) 

+ - 1)91/7 + 2fldqjl + 2-pdqP} {fldl + 2dq) 

= + 1) (dtfiijidi + 2dq) + dt]i{f^di + 2%)) . 

But {fidi + 2dq,X - iY) = and {yidi + 2dq, X - iY) = l/xp + 1, so that 
^'1^ = - n \2\,,2 ('^vV, X-iY). 

On the other hand 

dV _ _ 

— = ifxdi + 2dq)dtfi + (/iOi + 2dq)dtn , 

and we obtain 

dV _ {VvV,X-iY) {VvV,X + iY) iy^,,,m 

On replacing F with [/ = F/ 1 1 F 1 1 and the vectors X, Y with X = X/\\X\\,Y = 
Y/\\Y\\, respectively, we obtain 

^ = -^{VuU,X-iY){X + iY)-^{VuU,X + iY){X-iY) 

= -^fu. 

Conversely, if (7)(i) and (9)(ii) hold, then on reversing the above arguments, 
we obtain (7)(ii). 

In order to prove the last part of the Proposition, we calculate the time 
derivative of (9)(i), given that (9)(ii) holds. One of our aims is to see how 
the pair of equations (9) adapt to more general space times; this calculation 
is therefore carried out in Proposition 5.1 below, which establishes that in 
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the flat case, on writing a = {Cug){X — lY^X — iy), where now X,Y are 
supposed normahsed to have unit length, 

(ia{dt + U) = aa 

for some function a. In particular, by uniqueness of the solution of a first or- 
der PDE, if a vanishes on a hypersurface, it must vanish on the future causal 
space of this hypersurface. This completes the proof of the proposition. □ 



We now wish to see how, given a family of evolving unit vector fields 
Ut tangent to a conformal foliation, a correspond family of semi- 

conformal maps evolve. Since such a family of maps is only defined up to 
post-composition with a conformal transformation, we expect some gauge 
freedom in the choice of (ft- 

Corollary 3.3. Let Ut he a family of unit vector fields defined on evolving 
under the equation (9)(ii) and let ipt : Rf C be a corresponding family of 
(locally defined) semi- conformal mappings with fibres the integral curves of 
Ut- Then 

(10) d (^^) (U) = -t{^) . 

Furthermore, any other semi- conformal solution of equation (10) with the 
same fibres as (pt has the form fto^pt, where ft is an arbitrary conformal 
mapping of a domain of the complex plane C . 

Remark 3.4. The key point of the last part of the corollary, is that there 
is no restriction on the conformal mapping ft imposed by (10). 



Proof. Differentiation of the identity (iipt{Ut) = yields the equation 
By the fundamental equation (8) of a semi-conformal map, we have 

as required. 

Let ij^t be a family of semi-conformal mappings having the same fibres 
as ift so that locally, we can set ipt = ft° 'Pt for some family of conformal 
functions ft : A ^ C, where A is an open subset of C. We claim that 

(11) d (U) + T(V') = df (d (U) + T(^)] + AV(/) O ^ , 



p. BAIRD, M. WEHBE 



13 



where A denotes the dilation of 93. In particular, if (/? satisfies (10), then ^0 
satisfies (10) if and only if t(/) = 0; but this is so for ft arbitrary conformal. 

In order to show (11), we calculate d(^^^{U). First, 

Define <I> : / x ^ C (/ C R) by $(t,x) = (pt{x) and write for the 
connection in the pull-back bundle ip~^TC. Then d (^^^ {U) = and 
we obtain 

dt 



v^(d/(^)) = v^(d/od<i>(5/at)) 
= vX(d(/o$)(a/at)) 



u 

= vd(/ o $)([/, d/dt) + d(/ o <^){Vu{d/dt)) . 

The formula for the second fundamental form of the composition: Vd(/ o 
$) = d/(Vd^>) + Vd/(d$,d$) now gives 

Vfj (df (^)) = df (Vd^U, d/dt)) + Vdfid^U),d^d/dt)) + d(/ o <^){Vuid/dt)) 



d/(v*d$(a/at))=d/(v^^ 



d/(Vd$([/, d/dt) + d<i>{Vu{d/dt)) (since d<l>{U) = 0) 

dip 

'dt 

On the other hand, d o ip^ (U) = 0, so that 

Finally, the law for the tension field of a composition gives: 

r(0)=d/(r((^)) + AV(/), 
and (11) follows. □ 



We can normalise the evolution dpt/dt^ at least locally, in the following 
way. Let A C Rj] be an open subset on which pQ is submersive with con- 
nected fibres. Let S C A be a slice, that is, a 2-dimensional hypersurface 
which each fibre component intersects transversally in a unique point. For 
example, after a suitable change of coordinates, we could take S to be the 
slice: = 0, where we would then require dpQ/dx^{x) / for ah x G Sn^. 
Points of S would then be parametrized by the coordinates (x^,x^). 
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Let c(s) be an integral curve of U starting at a point of S: c(0) G 
S, c'(s) = U{c{s)). Then for sq sufficiently small that c(s) G A for all 
< s < So, we have 

(12) ^(c(.o)) = -£V^(c(.))ds. 

Note that r^(c(s)) G r^(c(s))C = Ty^C, where yo = ip{c{0)) = ip{c{s)) (since 
c(s) parametrizes the fibre of through c(0)), hence the integral on the 
RHS of (12) is well-defined and gives a vector in Ty^C = C. Of course, 

different choices of S will give different values for dipt/dt, but d (^^^^ {U) 
is well-defined and independent of this choice. 

Example 3.5. Let (po : — > C 6e the semi-conformal map ipo{x^ , x"^ , x^) = 
sj (x-'-)2 -|- (x2)2 + . Then the fibres of ipo are circles whose centres lie on 
the x^-axis lying in planes parallel to the x^x'^-plane. On letting Uq be the 
unit vector field tangent to the fibres of ipo, we have 

^ _ f x^ -x^ \ 

Take S to he the slice x"^ = 0; then an integral curve oJUq through {0,x'^,x^) 
is given by 

c(s) = (x"^ sin ( — ) , x"^ cos 



Now = l/^(xi)2 + (x2)2, so that 

In fact it is possible to calculate the complete evolution of ipo to be 

. o f r — it 
ipt = ix +r -targl — ^ 



where r = ^ {x^Y' + {x^Y' ~ '^^V is defined and smooth on the 

domain {x^)"^ + (^2)2 > t'^ , 



4 Coupled evolution of a metric and a unit vector field 

Consider a triple {M^, gt,Ut) consisting of a 3- manifold endowed with 
an evolving family of Riemannian metrics gt , together with an evolving unit 
vector field Ut- We suppose that the parameter t belongs to some connected 
interval / C R and we consider the space-time I x endowed with a metric 
g = -f{t,xfdt^ + gt, where / : I x ^ R is a non-vanishing smooth 
function. Then on / x M^, the vector field dt + fU is null; we will call 
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its integral curves the null curves generated by U. Our aim is to determine 
conditions when these are geodesic as well as either shear-free or generate 
conformal foliations. 

We first determine some basic quantities deriving from the metric G. In 
what follows, we fix a system of local coordinates (x" = t,x^,x'^,x^) and 
use the following range of indices: € {1, 2, 3}; a,b, . . . G {0, 1, 2, 3}. 

Recall that the Christoffel symbols are given by 



where we write R{da,db)dc = ^ di,dc - Va^Vg^^c = Ribc^d- We also use 
the convention G{R{da,db)dc,dci) = Robed = QdeR'ahc- ^^^^ ^^^o have oc- 
casion to consider the Christoffel symbols and curvature of the Riemannian 
metric gt] when there is the possibility of confusion, we will write ^T^^ or ^r.^^- 
to distinguish the two cases; similarly for the curvature and the connection. 

Lemma 4.1. With the notational conventions for indices given above, the 
Christoffel symbols of the metric Q = —f{t,x)'^dt'^ + gt are given as follows: 



Lemma 4.2. With the notational convention for indices given above and 
writing Vd/ for the second fundamental form of the function f considered as 
a function on (with t fixed), the components of the Riemannian curvature 
tensor of Q are given as follows: 



RoiOj = -\9^^dtgikdtgji + \dt{dtgij) - f\/df{di,dj) - \dt\iifdtgij 
^ijk = Tp {dij^t^ik - 9udt^]k - di{lnf)dtgjk + dj{ln f)dtgik^ 
^^ijk = ^^ijk + Tpidtgjkdtgu - dtgikdtgji) . 



Proof. In fact these follow from standard identities, the first known as the 
Mainardi equation, the second the Codazzi identity and the third, the Gauss 
equation, see [5]. We will prove the second of the above formulae, the other 
two being routine calculations. We have 



Kb = (daGbd + dbGad - ddGab) 



and the Riemannian curvature by 



rid a -pd r) _l F"^ "T*^ 

^abc ~ '^a.'- be '-'b^ ae ' be'- ae ae'- be ' 




= dtinf, r'], = aan/, r,o = fg'^d,f, 



R 



,0 

'ijk 




{di{hif)dtgjk 




dj {In f)dtgik) 



16 



SHEAR-FREE RAY CONGRUENCES 



However, note that 

diQjk - djQik = '^gij^lk + dkgij = -2guT'-jf. - dkQij 

so that 

diQjk - djQik = gij'^\k - 9li^]k ■ 
Substitution of this equahty into the above expression now gives the formula 
of the lemma. □ 



From Lemma 4.1, we have the following consequences. 

Corollary 4.3. The covariant derivative ofU and dt + fU along null curves 
defined by U is given respectively by the expressions: 

(i) "^i+fuU = ^ + f^fjU + lg'\dtguyd, 

+ (din f{U) + ^u'u^dtgi.j)dt 

(ii) Vl^fuidt + fU) = fgTadgf + {dt + fU)if)U 

+ ^u'g^''dtgijdk + (dt In / + Uif))dt . 

where, on writing U = u^{t,x^)di with respect to (fixed) local coordinates 
(x^) on M, we set dU/dt = {du'/dt)di . 



In what follows, we suppose that T = is a symmetric covariant 2- 
tensor defined on each {M^,gt) satisfying T(U,U) = 0, which, along with 
the function / is to be determined. In view of Corollary 4.3 (i), in order to 
preserve the length of U, we introduce the coupled evolution: 



(13) 



« f 



-f'^fjU-lg'''{dtgu)u'dk 
-2eQ{df\TM) + 2T 



where 9 = 6{t) is the 1-form dual to U: 9 = g{U, • ), a /? is the symmetric 
product of 1-forms: a Q f3 = ^{a /3 + /3 a) and U = u^di, where (x*) 
is a system of local coordinates on M^. Note that in the case of the flat 
metric, when we take / = 1 and T = 0, we recover the evolution (9)(ii). If 
we substitute into (i) the expression for dg/dt given by (ii), we obtain the 
alternative form: 



(14) 



(i) 
(ii) 



dU_ 

%_ 

dt 



-fVfjU+^gr^dgf + ^U{f)U 
-2 9Q{df\TM) + 2T 



where T is the linear transformation defined by g{T{E), F) 
each point {t,x) for all E,F £ T^M^. 



T{U) 



T{E,F) at 
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In general, the null vector field dt+fU is no longer geodesic, the expression 
for Vat+/f/(9j + fU) being given by Corollary 4.3, however, the following 
lemma gives the condition when it is geodesic. 

Lemma 4.4. Let the metric gt and the vector field Ut evolve according to 
equations (14). Then 

(i) {dt + fU)g{U, U) = 0, in particular, if U is of unit length on the 
hypersurface t = 0, then it remains so for all time; 

(ii) Vg^_|_jjj[/ = 0, in particular U is parallel along null geodesies deter- 
mined by U; 

(iii) {dt + fU)Q{U, dt) = 0, in particular, if U is tangent to for t = 0, 
then it remains so for all time; 



(iv) The null field dt + fU is geodesic if and only if 

^gradg/ + 



TiU) = --giadgf + -Uif)U, 



where T is the linear mapping determined at each point x by the identity 
g{f{E),F) = T{E, F) for all E,Fg T^M\ 

Proof. The first part of the Lemma follows directly from (13) (i) and the 
second from Corollary 4.3(i). To show part (iii), it suffices to note that 

= {dt\nf + U{f))dt + /gradg/ + ^-u' g''^ dtgtjdk , 

so that, from (ii), 

{dt + fu)g{u,dt) = g{u,vl^f^dt) = fu{f) + L^dtg){u,u) = o. 

For part (iv), the field dt + fU is geodesic if and only if 

Vd,+fu{dt + fU) = a{dt + fU), 
for some function a. By Corollary 4.3(ii), 

V|+^^(at + fU) = /grad,/ + {dt + fU){f)U + ^v'g^'dtgijdk + (dt In/ + Uif))dt 
= {dtlnf + U{f)}{dt + fU) + fgradgf 

+Lu'gkj^_0^d,f - ejdif + 2Tij)dk 

= {dt In / + Uif)} (dt + fU) + ^grad,/ - ^U{f)U + ff{U) . 
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Since the last three terms are tangent to (there is no ^t-component), the 
geodesic condition is equivalent to the vanishing of their sum, which gives 
the equation of (iii). □ 

Let (M^,(7j, Ut) evolve according to (14) and let {Xt,Yt} be an evolving 
orthonormal frame tangent to and orthogonal to Ut, which we call a 
complementary frame. Write Z = X — lY and let a = a{t, x) and p = p{t, x) 
{{t,x) G / X M^) be the quantities defined by: 

a = U^ug){Z, Z) = g(V|C/, Z) = Q{V%U, Z) 
p = g{V'-U,Z) = g{Vp,Z). 

The first of these is called the shear (see [14]) and measures the discrepancy 
from being conformal of the foliation determined hy U. It is quadratic in 
Z and so if we choose an alternative complex frame Z' = e'^Z, then the 
shear is now given by a' = e^'^cr, so that although a is not an invariant 
quantity, the property of its vanishing is. When necessary, we shall write 
cr{Z, Z) to indicate the dependence on Z. On the other hand, p is invariant 
independent of the choice of Z. Specifically, 

g(V|C/, Z) = -g{U, V^X + Vf.y - i[X, Y]) 

where g{U, V^X + Vf-F) is the mean curvature of the distribution on 
orthogonal to the vector field U, and where g{U, [X,Y]) measures the inte- 
grability of this distribution. 

We are mainly interested in the evolution of the shear a, however, we will 
find that both quantities a and p are intimately related and so we need to 
calculate both of their evolutions under the flow (14). Before doing this, 
we establish some basic formulae that will be useful in the computations 
that follow. In what follows, we assume that gt and Ut evolve according to 
(14). The following lemma is a simple calculation involving the Christoffel 
symbols as given by Lemma 4.1. 

Lemma 4.5. (i) For a vector field E = e^di tangent to M^, we have 
VEdt = E{lnf)dt - ^g{E, U)gTadglnf- ^E{lnf)U + e*r,,/^9fc ; 

g{VlZ, Z)=T{Z, Z); 

There are two concepts of shear-free ray congruence that we can consider 
on the space-time Ai^ = I x endowed with the metric G, both of which 
coincide in the flat case. In both cases, we suppose that the metric gt and 
the unit vector field Ut evolve according to (14) with the condition of Lemma 
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4.4(iv) satisfied, so the null curves generated by U are geodesic. The first 
notion requires that Lie transport of vectors in the 'screen space', that is 
the orthogonal complement of U in TM^, along the light rays be conformal; 
as in the flat case, we shall call this a shear-free ray congruence (SFR). 
The second concept requires that the evolution (14) preserves the property 
that the vector field U be tangent to a conformal foliation; we shall refer to 
the corresponding ray congruence as one that generates conformal foliations 
(CFGR). The following proposition characterises shear- free ray congruences. 

Proposition 4.6. Let {M^,gt, Ut) evolve according to (14) and suppose that 
T{U) = -^gradg/ + \U{f)U , so that the null curves generated by U are 
geodesic. Then the congruence is shear-free in = — /^di^ + gt) if and 

only if 

fa{Z,Z) + T{Z, Z) = 0, 
for any complementary frame Z . 



Proof. The condition that the ray congruence be shear-free is given by 

{^{dt+fu)Q) {z,z) = 0, 

equivalently by 

giVzidt + fU),Z) = 0. 
The condition of the proposition now follows from Lemma 4.5(i). □ 



In order to generate an SFR, a natural evolution now presents itself, given 
by the following theorem. Write /i = V^{7 for the mean curvature of the 
integral curves of U and let fi^ = g{fi, ■ ) be its dual. 

Theorem 4.7. Let f : I x ^ R (/ C R open, connected) be a non- 
vanishing function and let {AL^ ,gt,Ut) evolve according to the equations: 

( dU 

(i) 7^ = -fyfjU + gradgf 

\ da 

(ii) ^ = -2eQ{2df\TM-ffj^^-uif)e)-fCug. 

Then the null curves generated by U form a shear-free ray congruence. Fur- 
thermore, provided that is compact, for given initial conditions {go, Uq), 
the above system has a solution for small time. 



Proof. By the geodesic condition Lemma 4.4(iv) and the shear-free condi- 
tion given by Proposition 4.6, the tensor T in (14) is determined by the 
conditions: 

1 . „ 1 
T{Z,Z) = -fa{Z,Z) 



T{U) = --gradj + -[/(/)[/ 
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This implies 

T = -dfQ9 + u{f)e^ + feQfi^- ^£ug ■ 

Substitution into (14) now gives (15). 

To prove the existence of solutions, we note that, with respect to a sys- 
tem of local coordinates (x*) on M^, the system (15) is a first order sys- 
tem of differential equations in the coefficients and gij of U = u^di and 
g = gijdx^dx^ , respectively. Hence, on a compact domain, for given initial 
conditions, we are assured of a solution for small time by Cauchy's existence 
theorem, see for example [7] (10.4.5). □ 

Remark 4.8. In the above theorem, there is no constraint on the initial 
vector field Uq. At first sight, this appears to contradict the flat space case 
described in Section 3, where each vector field Ut is necessarily tangent to a 
conformal foliation (a = 0). However, in general the metric gt will evolve 
and produce a space-time with non-trivial curvature. In fact, if we take 
/ = 1, then 

^ = 20 /i^ - Cug 

where the right-hand side vanishes if and only if U is tangent to a Riemann- 
ian foliation. 



5 Congruences that generate conformal foliations 

We now address the second generalisation of an SFR on Minkowski space, 
that of a family of evolving conformal foliations. This turns out to be much 
more difficult, but has the advantage of providing a field equation somewhat 
akin to the Schrodinger equation on space-time. 

We suppose once more that we have a triple (M^, g^, Ut) of a 3-manifold, a 
metric and a unit vector field that evolve according to (14), where t belongs 
to some connected interval / C R. For convenience of notation, in what 
follows, we set r = g{VuZ, Z). The following proposition expresses how a 
and p evolve along null curves generated by U. 

Proposition 5.1. Let {M^ , gt,Ut) evolve according to (14), then the evolu- 
tions of a, p : / X C along null curves generated by U are given by 
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the formulae: 

{dt + fU)a = Qr(Z,Z) + /r-^(p + p))a-^T(Z,Z)(p-2C/(ln/)) 

+Z(ln /) - r(C/, Z)) - + T{U, Z)^ giVfjU, Z) 

(16) +Q{^R{du Z)U, Z) - fRicc\3{Z, Z) , 

{dt + fU)p = -^{fp + T{Z,Z))p-^{fa + T(Z,Z))a + jZilnf)Zilnf) 
(17) +^T{Z,Z)t - Ricci^(9i, U) - /Riccif (?/, U) . 

Proof. Unless stated otherwise, in what follows, the connection V is the 
Levi-Civita connection with repect to the metric G. Then 

{dt + fU)a = g{V9,+fuVzU,Z) + g{VzU,V9,+fuZ) 

= Q{^R{dt + fU, Z)U + VzVa,+fuU + V[9^^fu,z]U, Z) + G{VzU, Vg^+fuZ) 

= g {SR{dt + fu,z)u + V[9^+ju,z]U,z) +g{VzU,V9,+fuZ) 

since by Lemma 4.4(ii), V a^^fjjU = 0. We now evaluate the various quan- 
tities in this expression. 

The set of vector fields {Z,Z,U,dt/ f} form a basis of the complexified 
tangent space T'^{I x M'-^), with the properties g{Z, Z) = g(Z,^) = and 
g{Z,Z) = 2; the other vectors being mutually orthonormal. Thus we may 
write 

Va,Z = aiZ + [3{Z + ^iU + 5i{dt/f) , 
for some functions ai, 71, which we calculate as follows: 
ai = ^g{Va,Z,Z) = ^T{Z,Z) (Lemma 4.5(ii)), 

Pi = lg{V9,z,z) = o 

71 = g{V9,z, u) = -g{z,va,u) = g{z, fvfjU) 



Si = -g (Vd,z, = ^g {z,VaA) = jg{z,fg^'difdk) = g{z,gva.dgf) 

Thus 

Va,Z = \t{Z,Z)Z + fg{Z, VfjU)U + Z{f)% . 



Now set 

VuZ = a^Z + + 72C/ + h{dtlf) , 
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for functions 02 > /52 j 72 1 ) which are given by 

a2 = ^giVuZ,Z) = ^T, (32 = 
j2 = G{VuZ,U) = -G{Z,VuU) 

62 = -Q {VuZ, ^) = (Z, Vudt) = -iZ(ln /) + ir(C/, Z) (Lemma 4.5(i)) . 
Thus 

VuZ = \tZ - g{Z, VfjU)U - ^ (Z(/) - 2T{U, Z)) dt , 
and we obtain 

Va,+fuZ = \iT{Z, Z) + fT)Z + jj + 2TiU, Z)) dt . 

Proceeding in the same way, we obtain 

VzU = )^pZ +]^aZ-^ - 2T(U, Z)) dt , 

to give 

g (VzU, Va,+fuZ) = \{T{Z, Z) + fr)a + j - T{U, Zf^ . 

Similarly, we obtain 

[dt + fU, Z] = ^iT{Z,Z) +fT- fp)Z - ^faZ - \z{f)U - f{Z) , 

where we recall that T is the linear mapping determined at each point x by 
g{f{E),F) = T{E, F) for all E,F e T^M'-^. From this we deduce that 

'^[d,+fU,z]U = ^{T{Z,Z)+fT-fp)VzU-^faV^U-^Zif)VuU-Vf^^-^U, 
and so 

Q{'^[a,+fU,z]U,Z) = ^{fr-fip + p))a-^TiZ,Z)p 

-(^+r(f/,Z)) giVuU,Z). 

On the other hand 

Ricci^(Z,Z) = Ricci^'(X,X) - Ricci»(y,y) - 2iRicci ^(X, y) 

= gC^RiU, X)X, U) +g{<^R) {Y, X)X,Y) - g{3 R{U,Y)Y, U) 

-g{'R{X, Y)Y, X) - 2ig{9R{U, X)Y, U) 
= -g{3R{U,Z)U,Z). 
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On applying Lemma 4.2, we deduce that 

gi^R{U,Z)U,Z) = g{9R{U,Z)U,Z) + ^(dtg{Z,Uf -dtg{U,U)dtg{Z,Z)^ 

= -Ricd^iZ, Z) + j^(^(^^- T{U, Z)) ' + U{f)T{Z, Z) 

Combining the various expressions, now gives equation (16). The calculation 
of (dt + fU)p proceeds in a similar manner; we omit the details. □ 



The striking property about equation (16) is that the right-hand side is a 
quadratic form q = q{Z, Z) in the complex vector Z. Now any two quadratic 
forms on a l-(complex) dimensional space are proportional. In effect, if 
we replace Z by the vector field Z = e^^Z, then q{Z,Z) = e'^^^q{Z,Z) 
and similarly for any other quadratic form, so, since a = cr{Z, Z) is also a 
quadratic form, provided a{Z, Z) ^ 0, the equation 

q{Z,Z)=^a{Z, Z) 

is well-defined, independent of the choice of Z and determines a well-defined 
function 7 = j{t,x). However, we cannot make this assertion at points 
where a may vanish and so cannot claim in general that {dt + fU)a is 
proportional to o", which would allow us to conclude that if a vanished on a 
hypersurface, then it would vanish along integral curves of dt + fU . Instead 
we seek conditions when q vanishes. We can make this precise with the 
following consequence of Proposition 5.1 and Lemma 4.4. In what follows, 
we write TZ for the symmetric curvature tensor on {M^,gt) defined by 

7^ = -/Riccif + symQi^Ridt, -p,-), 

where {sym S){X,Y) = ^{S{X,Y) + S{Y,X)^ denotes the symmetrisation 
of an arbitrary covariant 2-tensor. Note that TZ depends on the ambient 
space-time manifold {M^,Q) and the vector field U. 

Theorem 5.2. Let {M^,gt,Ut) evolve according to the equations: 



(18) 



dU 

(i) — = -/Vf^CZ + gradg/ 



)t 
dg_ 

dt 



-2 0(d/|TA/) + 2r 



where t belongs to a connected interval I C H, T is a symmetric covariant 
2-tensor satisfying 

f(C/) = -lgrad,/ + ![/(/)[/ 

and where f : I x — > R is a smooth function. Then the integral curves 
of the vector field dt + fU are null geodesies with respect to the space-time 
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metric Q = —f'^dt^ + gt. If further 

(19) -^T{Z, Z) [p - 2C/(ln /)) + jZ{ln ff + n{Z, Z) = , 

where Z = X — iY is a complex vector field tangent to with X, Y, U an 
orthonormal triple, then ifUo is tangent to a conformal foliation on [M^,go), 
it remains tangent to a conformal foliation on {M^,gt) for all future time t. 



Note that although T(U, ■ ) is determined, we are free to choose T on the 
orthogonal complement of U in TM^. However, T is required to be a real 
tensor and the complex coefficient p = g^V-^U, Z) means we cannot solve 
(19) in its most general form. In what follows, we describe a number of 
special cases when solutions to (18) and (19) exist. 



Case 1: Stable points of the evolution. These are given by the equations: 

gradgln/ = VfjU 

T = OQdflTM 



(20) 



In particular, we must have U{f) = and on writing p, = 'VfjU for the 
mean curvature of the integral curves of U, necessarily ^ is a gradient and 
T = fO Q p^ ■ If we now make the additional assumption that Uq is tangent 
to a conformal foliation, then since gt and Ut are stable (independent of t), 
then Ut = Uq is also tangent to a conformal foliation and a must vanish 
identically. It then follows from (16), equivalently (19), that 

Ricci9(Z, Z) = Z{\nff . 

In fact, given a conformal foliation on a Riemannian 3- manifold, one can 
calculate the Ricci curvature in terms of invariants of the foliation, see [2], 
Proposition 6.1. One then obtains 

Ricci5(Z,Z) = ]^{C^g){Z,Z)- p\zf 

= V^^dln/(Z,Z)-Z(ln/)2. 

Combining these two expression for Ricci ^(Z, Z) we have a necessary con- 
dition on /: 

V^dln/(Z,Z) = 2Z(ln/)2. 

Conversely, given a vector field U = Uq which is tangent to a conformal 
foliation with mean curvature p = '^fjU = grad^/ a gradient. Then on 
defining T by (20), we have a stationary point of the evolution (18). To 
summarize: stationary points correspond to vector fields U whose associated 
mean curvature vector field p is a gradient; if U is initially tangent to a 
conformal foliation, then it remains so. 
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A special case is when U has integral curves geodesic. Then we can take 
/ = 1 and T = 0. In this case, if the integral curves of U also form a 
conformal foliation, then locally, U is tangent to the fibres of a harmonic 
morphism, that is a semi-conformal map which is also harmonic. It is to be 
noted that then, necessarily Ricci^(Z, Z) = 0, see [4]. Furthermore, if U is 
a vector field which is tangent to a conformal foliation, then the property 
that ^ be a gradient implies that there is a conformal deformation of the 
metric g which renders the integral curves geodesic (see [4], Lemma 4.6.6). 
The converse may not be true. 



Case 2: The case when (M^,(7o) has constant curvature. In this case, we 
take T = and / = 1. Then by (18), dg/dt = and gt = go for all 
t. Since OTyOt = 0, by Lemma 4.2, it follows that i?^.^ = 0, so that 
R{dt, Z,U, Z) = 0. On the other hand, by the hypothesis of constant 
curvature, Ricci^ = eg for some constant c and 

remains tangent to a conformal foliation. In the flat case, examples can be 
calculated from the twistor correspondence described in Section 3. 

Example 5.3. Consider the parametrized surface in CP^ given by 

icz 



{z,w) 



where c € C is non-zero. The incidence relation = ^xaa'V^ takes the 
form: 

—cz = vz — q 
-\\plw = —qz + u 
giving 



f + c 

The spinor field tangent to the corresponding SFR is given by 



l^A = P 



v + c 



for some complex-valued function p. A semi-conformal map tpt : C 
having the field U = ji^^{\z\'^ — l,2z) tangent to its fibres, is given by 

^t{xi,q) = ^^(1^1^ - 2xi(c + t) + (t + cf) . 

By Corollary 3.3, ipt is uniquely determined up to post- composition with a 
conformal mapping of the plane ft. If we take ft to be the dilation ft{w) = 
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^^w, then we obtain the semi-conformal mapping ipt = ft ° ft, with the 
same fibres of ipt : 

Mx) = - 2a^i(c + t) + {t + cf) . 

In the limit as t ^ oo, this converges pointwise to the semi-conformal map- 
ping = I having geodesic fibres. This shows that in the limit as t ^ oo, 
the corresponding foliation given by the integral curves of U converges to a 
conformal foliation by geodesies. This is as we would expect when conver- 
gence as t ^ CO occurs, since such a foliation is stable with respect to the 
evolution (18). 

Example 5.4. For non-zero integers k and I, there is a smooth semi- 
conformal map if : ^ S"^ between Euclidean .spheres of the form 

¥p(cos(s)e^",sin(s)e''') = (cos a(s), sin a(s)e'('='^+'^) , 

where a : [0,7r/2] — > [0, tt] satisfies a(0) = 0, a{7r/2) = vr [4], Example 
13.5.3. If we use coordinates (x^,x^,x^) = {s,a,b) on 5^, then the canonical 
metric has the form 

g = ds^ + cos^ s da^ + sin^ s db^ . 

With respect to the basis d/ds, d/da, d/db of the tangent space, the non-zero 
Christoffel synbols are given by 

= sin s cos s, = — sin s cos s, = ~ tan s, = cot s , 

and the unit vector field tangent to the fibres of (p is given by 

Uo = ^ = {-Ida + kdt) . 

V k^ sin s + cos^ s 

Let us suppose that for each t, the vector field Ut has the form Ut = uds + 
vda + wdb for functions u, v, w which depend on s and t only. Then a cal- 

dU 

culation shows that the equation = —V^Ut is equivalent to the system: 

-— = —iv—w smscoss — n— — 
\ 

ov ov 
— = — u— — |-2u'ytans 

ow ow 
— — = —u— zuw cot s , 

Ot OS 

with initial conditions: 

u{s,0)=0, ^(s,0) = -=_L==, w{s,0) = ^==L=== 

V sm s + cos^ s V sm s + cos^ s 

The theory of partial differential equations of first order guarantees the ex- 
istence of a solution for small time. Note that, even though Uq is tangent to 
the fibres of the map (p : ^ S'^, for later time, we cannot expect that there 
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will be a globally defined semi-conformal mapping on with Ut tangent to 
its fibres. We can only be sure that locally such a mapping exists. 



Case 3: The case of integrable complementary distribution. If, via the flow 
(18), we evolve through vector fields tangent to a conformal foliation, then 
remarkably, the integrability of the complementary distribution is preserved, 
as described by the following lemma. 

Lemma 5.5. Let {M^,gt,Ut) evolve according to (18), then 

{dt + fU){p -p) = -^TiZ,Z)ip -p)- ^T(Z,Z)a + ^T{Z, Z)a . 

In particular, if p — 'p = 2ig{U, [X,Y]) vanishes at t = and a = 0, then 
p — J) = for future time. 



Proof. The formula is a direct consequence of (17), under the assumption 
of Theorem 5.2, that T{U) = — ^gradg/ + ^U{f)U. The conclusion then 
follows as a consequence of the uniqueness of solutions of a first order ODE 
for given initial conditions. □ 



On taking / = 1, this enables us to give a satisfactory form of the evolu- 
tion (18) in the case when initially the complementary distribution is inte- 
grable. 

Theorem 5.6. Let [AI"^, gt,Ut) evolve according to the equations: 
dU 



(21) 




,t -^'-^ 

J = -^{n-2n{u,-)Qe + n{u,u)e^] 



where v = g{U, V^X + Vf.y) is the mean curvature of the horizonal dis- 
tribution (C/"*") which we assume non-zero at all points of M^, and where 
TZ = — Ricci^ -\- sym.Q(^ R{dt, ■ )U, ■ ) . Then ifUo is tangent to a conformal 
foliation on {M^,go) with integrable complementary distribution U'^ , it re- 
mains tangent to a conformal foliation on (M^,gt) for all future time t for 
which p remains everywhere non-zero. Furthermore, the complementary dis- 
tribution remains integrable and the integral curves of the vector field dt + fU 
are null geodesies with respect to the space-time metric Q = —f'^dt^ + Qt- 



Proof. The result is a consequence of Theorem 5.2, on taking / = 1 and T 
to be the tensor 

-- (Riccif + symg(^ii:(5t, ■)U, ■)) =--7^, 
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normalised in order that T(U) = 0. For, with these assumptions, the require- 
ment (19) is satisfied. However, we require that T remain real, equivalently, 
that p remain real. Now by (16), Lemma 5.5, and our choice of / and T, we 
have the coupled system: 

r {dt + fU)a = {lT{Z,Z)+T-^{p + p))_a_ 

\ {dt + fU){p-p) = -lT{Z,Z){p-p)-^T{Z,Z)a + ^T{Z,Z)a, 

with initial conditions a = and p — J) = 0. Once more, by the uniqueness 
of solutions of a first order ODE for given initial conditions, it must be the 
case that a and p — J) remain zero for future time. We now set i/ equal to 
the real quantity —p. □ 

Remark 5.7. We can deal with the case when p = in a similar way, but 
now choosing f to be a function such that U{f) is nowhere vanishing (where 
this is possible). Then equation (19) requires that we take T to be the tensor 

normalized so that T{U) = 0. 

It is interesting that (21) (ii) is a modified Ricci flow (depending on the 
influence of the term Q{^R{dt, ■)U, •) in TZ), coupled to the evolution of 
the vector field U. Some examples of Ricci solitons evolve under the Ricci 
flow through semi-conformal mappings [2]. To prove existence of solutions 
to the system (21), even for short time, would seem to be quite tricky. 
The term Q{^R{dt, ■ )U, ■ ) that appears in the tensor TZ in general contains 
terms involving the derivative of the metric g; this will affect the nature 
of the equations. The following example gives an indication of how this 
might happen. Note that the example, consisting of a conformal foliation 
by geodesies, also fits into Case 1, and with the choice / = 1, T = 0, would 
correspond to a stationary evolution. However, with the imposition of the 
evolution of the metric given by Theorem 5.6, the metric will evolve in time. 

Example 5.8. We consider a radial metric on a domain o/R^ of the form 

g = a{t, rfdr"^ + b{t, rfgs2 = a^dr^ + b'^{ds'^ + sin^ s dO'^) , 

where the coordinates correspond via {x^,x'^,x^) = {r,s,9), with {s,9) stan- 
dard coordinates on S*^. Then the r -curves form a conformal foliation by 
goedesics with unit tangent given by U = The non-zero Christoffel 

symbols are given by 

T\, = 9,lna, = T^^ = -Mb!^, = dAnb, Ffg = dAnb, 

rig = — sin s cos s, rL = cot s . 
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The tensor TZ then has the following non-zero components: 
TZu = 2 (5„.^ In 6- a^. In a Sr. In 6 + (5r In 6)2) 
'^13 = cot s dr In b 

7^22 = -l + -dr(^] +-(dtr^b-dtlnadrb) 

a \ a J a ^ ' 

7^33 = sin^ s 7^22 • 

An orthonormal triple containing U is given by {X,Y,U} = {^^; fesLa M' ^J"' 
A calculation now gives 

V = g{U, V'^X + V^Y) = - ^drb . 

The evolution equation (21) is then just the single equation: 

'bdr-b 
a 

Stationary points, i.e. time-independant solutions require 

But now we can suppose that r is parametrised so that a{r) = 1. Then up 
to replacing r by r + c for a constant c, we have b = r, which corresponds 
to the Euclidean metric. 



bdfJ^b — dfb drb — bdt In a drb = a — dr 



6 The corresponding field equations and conformal invari- 
ance 

We conclude with some general remarks concerning the application of the 
above to the description of massless fields. Our perspective has been to 
couple the evolving field, which we suppose represented by the vector field 
U (or the corresponding semi-conformal map f), with the gravitational field 
g. Thus the two should evolve simultaneously following some well-defined 
equations. In order to arrive at these equations, we have adapted the flat 
case (Minkowski space) to more general curved metrics. Prom the physical 
point of view, we do not wish to "presuppose" the space-time manifold 
{Ai^,G), but rather to introduce a natural dynamic into the pair {gt,Ut) 
based on the wish to preserve the conformality of the foliation and the 
geodesity of the null curves determined by U. As such, we would like to be 
able to see semi-conformal mappings as deriving from a variational principle. 
In fact, this is the case if the mapping has values in a surface. 

Let if : (M"*, g) {N^, h) (m > 2) be a smooth mapping of a Riemannian 
manifold into a Riemannian surface. Then there are at most two non-zero 
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eigenvalues Ai^, of the first fundamental form ip*h, which we can suppose 
satisfy Ai > A2 > 0. Then 99 is semi-conformal if and only if Ai = A2 at each 
point. On the other hand, the energy density is the function: 



where dip* denotes the adjoint of the derivative, and the Jacobian the func- 
tion: 



Clearly e^p — j,p > with equality if and only if Ai = A2, i.e. if and only if cp 
is semi-conformal. If M™" is now compact (otherwise we restrict to compact 
domains), it follows that (p is semi-conformal if and only if it is absolutely 
minimizing with value zero for the functional: 



where dvg is the volume element on M™ and where the power m/2 assures 
the conformal invariance of / with respect to conformal deformations of g. It 
is not known whether in any given homotopy class, the functional / admits 
an absolute minimum and if so whether it is always zero. 

We can obtain a time-dependent differential equation governing the evolu- 
tion of in a similar way to Corollary 3.3 for the flat case. On differentiating 
the equation dipt{Ut) with respect to t and applying (18)(i), we obtain 



the latter equality coming from the fundamental equation (8) for a semi- 
conformal map. We can once more integrate this along integral curves of U 
to obtain: 



where the start point c(0) lies on a chosen hypersurface of which the 
integral curves of U hit transversally and x = c{so). We can view this equa- 
tion as a modified Schrodinger equation, or heat equation, with potential 
given by the function /. 

A further observation is as follows. If U is tangent to a conformal foliation 
on a 3-manifold {M^,g), now possibly with singularities, and ip : — > N'^ 
is a semi-conformal map with U tangent to its fibres, then where it is smooth, 
the vector field X'^U is divergence free, where A denotes the dilation of ip. 




= \/det(dv9 o dtp*) = A1A2 . 





dp^tiUt) = diptiVfjU) - diptigiad gf) 
T{ipt) - diptigiadgf) ; 
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In particular, if 5" C is a closed surface which avoids singularities of U 
(and of (f), then by Stokes Theorem, the quantity: 

Q{S) = [ X^g{U,n)dvs, 
Js 

is independent of the homology class of S in \ S, where S represents the 
singular set of U [17]. Furthermore, it is easily checked that the integral 
Q{S) is conformally invariant with respect to conformal deformations of 
the metric g. For example, if ip : — > is radial projection given by 
X 1-^ with a singularity at the origin and S is taken to be any sphere 

surrounding the origin, then Q{S) = An. We then expect Q{S) to represent 
some intrinsic quantity associated to the field given by if. 

One of the features of a shear-free ray congruence on Minkowsi space is 
its conformal invariance. Thus if ^u"^ satisfies equation (2) with respect to 
the Minkowski metric ^m*; then it does so with respect to any conformally 
equivalent metric 0,'^g^4, where : — > R is a smooth non- vanishing 
function. Our equations (15), (18) and (21) do not have this general invari- 
ance property. This is because of our requirement that the vector field U be 
parallel along null curves generated by U, which is not a conformally invari- 
ant property. Explicity, if we consider a conformal deformation Q = Vl^Q of 
the space-time metric Q = —f'^dt^ + gt, then for arbitrary vector fields, we 
have the relation: 

v|f = V|F + E{ln n)F + F(ln n)E - Q{E, F)grad g In O . 
If we now consider the new quantities U = U/Q, f = ^If, then 

"^l+fuU = ^^i+fuU + ^U{lnn)idt + fU) - ^gradg lnl7 , 

so that both U and U are parallel along the null curves they generate (with 
respect to Q and Q, respectively) if and only if 

r (dt + fum = 

\ U{n)U -gradgQ = 0. 

The second equation is quite restrictive. It says that the gradient of 
must be parallel to U, but then if Q is non-constant, this means that the 
complementary distribution JJ-^ is integrable. In particular, if [/-*- is not 
integrable, then the only conformal deformation possible has il. constant. 
On the other hand, if [/-*- is integrable, then any smooth nowhere vanishing 
function O which is constant along the integral surfaces of the distribution 
IJ-^ which also satisfies {dt + fU){Q) = (O is constant along the null curves 
generated by U) is an allowable conformal deformation. It is easily checked 
that the equations (15), (18) and (21) are invariant with respect to such 
conformal deformations. 
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However, we take the view that such a global conformal invariance should 
not be a requh'ement of our model. The concept of conformal foliation, 
or semi- conformal map, is a conformal invariant, so we can modify our 
initial conditions (M^, ^tq, C/q) by any conformal deformation: {M^,gQ = 
lo^Qq, Uq = ^Uq) and then allow the approprate evolution to take place. This 
is consistent with our perspective that the triple {M^,gQ, Uq) has an implicit 
dynamic associated to it. We could then view the conformal invariance of 
the initial conditions as a local conformal invariance. 
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